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ABSTRACT 



We have investigated the possibility of direct tests of little Higgs models incorporating triplet 
Higgs neutrino mass mechanism at LHC experiments. We have performed Monte Carlo studies 
of Drell-Yan pair production of doubly charged Higgs boson followed by its leptonic decays 
whose branching ratios are fixed from the neutrino oscillation data. We propose appropriate 
selection rules for the four-lepton signal, including reconstructed taus, which are optimized for 
the discovery of with the lowest LHC luminosity. As the Standard Model background can 
be effectively eliminated, an important aspect of our study is the correct statistical treatment 
of the LHC discovery potential. Adding detection efficiencies and measurement errors to the 
Monte Carlo analyses, <l>^^ can be discovered up to the mass 250 GeV in the first year of LHC, 
and 700 GeV mass is reachable for the integrated luminosity L = 30 fb~^. 
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1 Introduction 



The main motivation of the Large Hadron Colhder (LHC) experiment is to reveal the secrets 
of electroweak symmetry breaking. If the hght standard model (SM) Higgs boson H will be 
discovered, the question arises what stabilizes its mass against the Planck scale quadratically 
divergent radiative corrections. The canonical answer to this question is supersymmetry, pre- 
dicting a very rich phenomenology of sparticles in the future collider experiments. 

Alternatively, the light SM Higgs boson may signal some strong dynamics at high scale A ~ 47r/, 
where / is the decay constant of the new strongly interacting theory [1]. The most interesting 
class of models in such a scheme are the little Higgs models [21 El H]. In those models the SM 
Higgs boson is a pseudo Goldstone mode of a broken global symmetry and remains much lighter 
than the other modes of the model, thus solving the little hierarchy problem and postponing 
the solution to the fundamental hierarchy problem to the scale A. Those models are also very 
interesting from collider physics point of view since they predict the existence of new particles, 
such as a new set of heavy gauge bosons Wh, Zh, a vectorlike heavy quark pair T, T with 
charge 2/3, and triplet Higgs bosons If the new particle masses are 0{1) TeV, direct tests 
of the models are possible at LHC [51 El [7] . 

An important open issue to address in the context of little Higgs models is the origin of non- 
zero neutrino masses [HI El [IDl El [12]. The neutrino mass mechanism which naturally occurs 
in those models is the triplet Higgs mechanism [131 113] (sometimes called type II seesaw) which 
employs a scalar with the SU (2)^ x f/(l)y quantum numbers $ ~ (3, 2). The existence of such a 
multiplet in some versions of the little Higgs models is a direct consequence of global symmetry 
breaking which makes the SM Higgs light. For example, in the minimal littlest Higgs model [15] . 
the triplet Higgs with non-zero hypercharge arises from the breaking of global S'f/(5) down to 
SO(b) symmetry as one of the Goldstone bosons. Its mass M$ ~ ggf, where Qs < 47r is a model 
dependent coupling constant in the weak coupling regime P], is therefore predicted to be below 
the cut-off scale A, and could be within the mass reach of LHC. Although the triplet mass scale 
is 0{1) TeV, the observed neutrino masses can be obtained naturally. Firstly, non-observation 
of rare decays /i eee, n — ^ 67, r — > where £ = e, /i, implies that the triplet Higgs boson 
Yukawa couplings Yij must be small, thus suppressing also the neutrino masses. Secondly, the 
vaccuum expectation value (vev) of the neutral component of triplet vq> contributes at tree level 
to the SM oblique corrections, and is therefore severely constrained by precision data. There 
exist additional mechanisms which can explain the smallness of f $ in little Higgs models. Since 
the smallness of w<i> is the most natural explanation of the smallness of neutrino masses in the 
little Higgs models, we assume this to be the case in this work. 

The aim of this paper is to study the possibility of direct tests of little Higgs models and 
neutrino mass mechanisms at LHC experiments via pair productions and subsequent decays of 
triplet Higgs boson. We study the Drell-Yan pair production of doubly charged component of 
the triplet [ISl [IH dHl [ISl[2ni[2I] 

pp $++$ , (1) 

followed by the leptonic decays. Notice that (i) the production cross section does not depend 
on any unknown model parameter but the mass of {ii) smallness of in this scenario, 

following from the smallness of neutrino masses, implies that the decays — > are 
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negligible, and we neglect this channel in the following analyses; (iii) the leptonic decay 
branching fractions do not depend on the size of the Yukawa couplings but only on their ratios 
which are known from neutrino oscillation experiments. In the triplet model the normally hier- 
archical light neutrino masses predict BR{^^^ — > n'^) ~ BR{^^^ — > r+r+) ^ — >• 
/i+r+) 1/3. Therefore this scenario is predictive and testable at LHC experiments. 

The production process ([T]) has been studied before in various theory papers. In this work we 
first carry out a pure Monte Carlo study of the signal and background processes in the environ- 
ment of LHC detectors. After that we improve our analyses by adding particle reconstruction 
efficiencies and Gaussian distortion functions for particle momentas and EJp^''^'^. Those mimic 
the detector inefficiency effects at the Monte Carlo level. We believe that those effects help us 
to estimate the realistic mass reach of the LHC detectors to the process under study. 

In our study the new results are the following. For the signal reconstruction we use new criteria, 
such as equality of invariant masses of positively and negatively charged leptons together with 
total UpT cut for all leptons, which allows us to achieve better reconstructions efficiencies 
compared to the standard cuts. We also reconstruct tau lepton final states with more than one 
r, which has not done before in this context. As all the SM background can be eliminated in the 
case of this process, correct statistical analyses of the results in the limit of no background is an 
important aspect of our study. For the discovery criteria we have used the Log-Likelihood Ratio 
(LLR) statistical method to demand 5a discovery potential to be bigger than 95% (1 — CL^+j, > 
0.95). Our results are optimized for the discovery of process ([T]) with the lowest possible LHC 
luminosity. The pure Monte Carlo study shows that ^'^'^ up to the mass 300 GeV is reachable 
in the first year of LHC (L = 1 fb~^) and up to the mass 800 GeV is reachable for the 
luminosity L = 30 fb~^. Including the Gaussian measurement errors in the Monte Carlo the 
corresponding mass reaches become 250 GeV and 700 GeV, respectively. The errors of our 
estimates of the required luminosity for discovery depend strongly on the size of statistical 
Monte Carlo sample of the background processes. 

The paper is organized as follows. In Section 2 we present the collider phenomenology of triplet 
Higgs boson and relate collider observables to neutrino mass measurements. In Section 3 we 
discuss the Monte Carlo produced signal and background processes. In Section 4 we present 
the details of reconstruction and analysis procedure and results. Detector effects are discussed 
in Section 5. Finally we conclude in Section 6. 



2 Neutrino masses and collider phenomenology 

In this work we consider little Higgs scenarios in which, due to the breaking of global symme- 
try protecting the SM Higgs boson mass, the spectrum of the model contains also a pseudo 
Goldstone boson with the SU{2)l x U{1)y quantum numbers $ ~ (3,2) [ISII22]- Although $ 
is predicted to be heavier than the SM Higgs boson, the little Higgs philosophy implies that its 
mass could be 0{1) TeV [1]. Due to the specific quantum numbers the triplet Higgs boson cou- 
ples only to the left-chiral lepton doublets Li ~ (2, —1), i = e, fj,, r, via the Yukawa interactions 
given by 

L = iL'iT2Y'^{T-<^)Lj + h.c., (2) 
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where Yij are the Majorana Yukawa couphngs. The interactions (E]) induce lepton flavour 
violating decays of charged leptons which have not been observed. The most stringent constraint 
on the Yukawa couplings comes from the upper limit on the tree-level decay fi eee and i^ 
YeeYefj, < 3 ■ 10~^(M/TeV)^ [m EH]. Experimental bounds on the tau Yukawa couplings are 
much less stringent. In our collider studies we take Yrr = 0.01 and rescale other Yukawa 
couplings accordingly. In particular, hierarchical light neutrino masses imply Yge, Y^^ ^ VVr 
consistently with the direct experimental bounds. 

According to Eq.(l2]), the neutral component of the triplet Higgs boson $° couples to the left- 
handed neutrinos with the same strength as couples to the charged leptons. If $° acquires 
a vev non-zero Majorana masses are generated for the left-handed neutrinos [131 E]. Non- 
zero neutrino masses and mixing is presently the only experimentally verified signal of new 
physics beyond the SM. In the triplet neutrino mass mechanism the neutrino masses are given 
by 

{m^)ij = YijVq,. (3) 

We assume that the smallness of neutrino masses is explained by the smallness of f In a realis- 
tic scenario massless Majoron, the Goldstone boson of broken lepton number, must be avoided. 
This is achieved by an explicit coupling of $ to the SM Higgs doublet H via ^^^H^H^ [H], 
where n has a dimension of mass. If ~ M$, in the concept of seesaw [26] the smallness of 
neutrino masses is attributed to the very high scale of triplet mass M$ because f$ = fiv^/M^, 
where v = 174 GeV. However, in the little Higgs models the triplet mass scale 0{1) TeV alone 
cannot suppress Therefore in this model fi <^ M$, which can be achieved, for example, via 
shining of explicit lepton number violation from extra dimensions as shown in ref . fI7[ EH] , or 
if the triplet is related to the Dark Energy of the Universe [291 EO]. Models with additional 
(approximate) T-parity [22] make the smallness of f $ technically natural. However, if the T- 
parity is exact, f$ must vanish. In this work we do not consider the naturalness criteria and 
assume that the above described neutrino mass scenario is realized in nature. In that case 
Yvi^, ~ 0(0.1) eV while the Yukawa couplings Y can be on the order of charged lepton Yukawa 
couplings of the SM. As a result, the branching ratio of the decay $ — > WW is negligible. 
We also remind that f $ contributes to the SM oblique corrections, and the precision data fit 
T < 2 • 10~^ [31] sets an upper bound f$ < 1.2 GeV on that parameter. 

Notice the particularly simple connection between the flavour structure of light neutrinos and 
the Yukawa couplings of the triplet via Eq.Q. Therefore, independently of the overall size 
of the Yukawa couplings, one can predict the leptonic branching ratios of the triplet from 
neutrino oscillations. For the normally hierarchical light neutrino masses neutrino data implies 
negligible $ branching fractions to electrons and BR{^~^^ — »• /x^/i^) ~ i?i?($++ — > t~^t^) ^ 
BR{^~^~^ — i> /x^r+) ^ 1/3. Those are the final state signatures predicted by the triplet neutrino 
mass mechanism for collider experiments. 

At LHC can be produced singly and in pairs. The cross section of the single production 
via the WW fusion process [H] qq scales as ~ f|. In the context of the littlest Higgs 

model this process, followed by the decays W~^W~^, was studied in ref. [3 [71 [32]. The 

detailed ATLAS simulation of this channel shows [32] that in order to observe an 1 TeV 

-•^In little Higgs models with T-parity there exist additional sources of flavour violation from the mirror 
fermion sector . 



3 



one must have > 29 GeV. This is in conflict with the precision physics bound w$ < 1.2 GeV 
as well as with the neutrino data. Therefore the WW fusion channel is not experimentally 
promising for the discovery of doubly charged Higgs. 

On the other hand, the Drell-Yan pair production process pp is not suppressed 

by any small coupling and its cross section is known up to next to leading order |19j (possible 
additional contributions from new physics such as are strongly suppressed and we neglect 
those effects here). Followed by the lepton number violating decays i^i^, this process 

allows to reconstruct invariant mass from the same charged leptons rendering the SM 
background to be very small in the signal region. If one also assumes, as we do in this work, 
that neutrino masses come from the triplet Higgs interactions, one fixes the leptonic 
branching ratios. This allows to test the triplet neutrino mass model at LHC 

3 Monte Carlo simulation of the signal and backgrounds 

The production of the doubly- charged Higgs is implemented in the PYTHIA Monte Carlo 
generator [33]. The final and initial state interactions and hadronization have been taken into 
account. We have used the CTEQ5L parton distribution functions. 

In the following analysis the normal hierarchy of neutrino masses and a very small value of the 
lowest neutrino mass is assumed. Such a model predicts that Higgs decay into electrons can be 
neglected and that there are three dominant decay channels for with approximately equal 
branching ratios: 



• 


$±± - 




• 


$±± - 




• 


$±± - 





We have studied only pair production of doubly charged Higgs due to the reasons pointed out 
above. pair decay products can combine to five different r and /i combinations: 4/i, 3/ilr, 
2/x2r, l/x3r and 4r. Before reaching the detector, r decays into an e, /x or a hadronic jet 
(marked as j below) with branching ratios of 0.18, 0.17 and 0.65, respectively [3l]. r hadronic 
jets and /i-s are well visible and reconstructible in detector. The reconstruction of an energetic r 
from electron decay is sensitive to detector effects, involving sophisticated background processes 
|35j . In the current analyses we will neglect this channel, which will cause 31% loss of the total 
signal. Such loss is still sufficiently low and can be considered acceptable. Table [1] gives the 
cross sections and the Monte Carlo generated event numbers in our study. 

The signatures of $ decay are very clean due to (i) high transfer momentum of the decay 
products, [ii) lepton number violation and {iii) pair production of $. The Standard Model 
particles are lighter than $, so the background /i-s and r-s must have smaller transverse energy 
and they do not produce an invariant mass peak in /i'*'/i^, yU^T+, r+r+, final states. The present 
lower bound for the invariant mass of $ is set by Tevatron to M$ > 136 GeV [3S| [37] . In our 
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Process 


Total (T 
ifb) 


N of events 
generated 


Corresponding 
luminosity (fb^^) 


Signal 

M$=200 GeV 
M$=500 GeV 
M$=1000 GeV 


7.78E+01 
1.99E+00 
5.58E-02 


l.OOE+05 
l.OOE+05 
l.OOE+05 


1.28E+03 
5.03E+04 
1.79E+06 


Background 
pp^tt ^ U 
pp ^ ti Z 
pp — * ZZ 


8.84E+04 
6.50E+02 
2.12E+02 


2.55E+07 
1.50E+05 
l.OOE+05 


2.88E+02 

2.3E+02 

4.72E+02 



Table 1: Cross-sections, numbers of Monte Carlo generated events and the corresponding inte- 
grated luminosities of the generated events. For the signal events we have taken the branching 
ratios BR($±± ^ = BR($±± ^ /i^r±) = BR(<I>±± ^ r±r±) = 1/3. 

study, four-lepton background processes with reasonable cross-sections and high px leptons 
arise from three Standard Model processes 

• pp ^ ti, 

• pp ^ tiZ, 

• pp ZZ . 

PYTHIA was used to generate ti and ZZ background {tt is forced to decay to WWbb and W 
leptonically). The CTEQ5L parton distribution functions were used. CompHEP was used to 
generate the Ztt background via its PYTHIA interface [3Hl[M]- All the datasets were generated 
in Baltic Grid. In addition to background processes shown in Table [H some other four-lepton 
background processes exist involving 6-quarks in the final state (for example, pp —>■ bb). As 
such processes are very soft, it is possible to use the effective tagging methods [IQ] and totally 
eliminate this soft background [H]. Also, we do not consider possible background processes 
from the physics beyond the Standard Model. 

4 Reconstruction and analysis of the Monte Carlo data 

To study the feasibility of detecting the signal over background, we have to work with five 
possible reconstruction channels according to the following final states. 

• 4/i: The cleanest and most simple channel. 

• 3yulr: The channel is easily reconstructable using an assumption that the 
neutrino originating from the r decay is coUinear with r-jet and gives majority to the 
missing transverse energy (i?^*'^*). 
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• 2/i2r: The signature can be reconstructed using the same assumptions for 
both r-neutrinos. The whole EJp^'^'^ vector has to be used here, while in the previous 
channel only one component was needed. 

• l|u3r: The channel can be reconstructed theoretically relying on an addi- 
tional requirement that the two Higgs bosons have equal invariant masses. However, the 
reconstruction is very sensitive to the experimental accuracy of i?™*'^ determination. 

• At: The channel can not be reconstructed (and triggered by the single muon 
trigger) . 

First, we apply general detector related cut-offs for the Monte Carlo generated data. Generated 
particles were reconstructed within the pseudorapidity region |?7| < 2.4 and with transverse 
momentum higher than 5 GeV. These are the natural restrictions of the CMS and ATLAS 
detectors at the LHC. Only the pseudorapidity region \ri\ < 2.4 is reachable for the detector 
and only the events with > 5 GeV are typically triggered. These restrictions suppress 
mainly the soft Standard Model background. The efficiency of lepton reconstruction and charge 
identification rate are very high, we use the values 0.9 and 0.95, respectively [12]. 

The invariant mass of two like-sign /i-s and/or r-s are calculated using equation: 

imf^r = {pt+ptr, (4) 

where pi^2 is the /x or r 4-momentum. Since the like-sign signal of fi-s or r-s originate from 
a doubly charged Higgs boson, the invariant mass peak measures the mass of doubly charged 
Higgs, mj = M$. 4-/i final state allows to obtain invariant masses directly from Eq.(jl]). In 
channels involving one or several r-s, which are registered as r-jets or secondary /x-s (marked 
as fi' below), the momenta of jets has to be corrected according to the equation system: 



Pr = k^I^et, (5) 

PTmiss = ^PTi^) (6) 

i 

M$++ = M$— , (7) 



where i counts r-s, p marks 3-momentum, p^^ is the vector of transverse momentum of the 
produced neutrinos, PTmiss is the vector of missing transverse momentum (measured by the 
detector) and ki > 1 are positive constants. Eq.(l5]) describes the standard approximation 
that the the decay products of a heavily boosted r are collinear |12]. Eq.(l6l) assumes missing 
transverse energy only to be comprised of neutrinos from r decays. In general, it is not a 
high-handed simplification, because the other neutrinos in the event are much less energetic 
and the detector error of E^'''^^ is order of magnitude smaller [l3]. Using the first two formulas, 
it is possible to reconstruct up to two r-s per event. Additional requirement of Eq.(I7]) allows 
to reconstruct the third r per pair event, although very low measurement errors are needed. 

A significant fraction of r-s (0.18) decay into yu'-s that cannot be distinguished from primary 
lis in the detector. Still, if reconstructed invariant masses of and $ are considerably 
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different, we can suspect tliat one or several fi-s originate form r decays. In sucli case we can 
again use Eq.([S])-([7j) to correct tlie 4- momenta of decay products. Wlien only one secondary 
muon is present, i?™'''' points into tlie same direction as its px- Otherwise E^^^''^ is a superpo- 
sition of neutrino transverse momenta. Such correction tightens the invariant mass peak of the 
signal and does not produce any artificial background. 

The occurrence probability of different reconstruction channels are presented in Table [21 The 
second column shows probabilities of Higgs decay to /i-s and M r-s. Next columns describe 
the final state after r decay to /x'-s and/or jets. Different columns mark the number of secondary 
fi-s and the rows designate r-jets in the detector recordings. The events having at least one 
r — > e in a final state are omitted in our analysis, as well as events with M > 3 or iV(/i') > 2. The 
proportions of reconstructible signatures are marked in a bold-face. The table shows that 0-3 
jet channels together with /i correction are almost equally important and overall reconstructible 
channels comprise 64% of total events. 



Decay channel 


After Higgs 




After 


T decay ( x = 


= jetr) 






decay (x = r) 


(r ^ /i') 


1 (r ^ /i') 


2 (r ^ /i') 


3 (r /i') 


4 (r /i') 


2$ ^ 4/x 


0.1111 


0.1111 


0.0377 


0.0107 


0.0012 


0.0001 


2$ 3/ilx 


0.2222 


0.1443 


0.0736 


0.0125 


0.0014 




2$ 2/i2x 


0.3333 


0.1407 


0.0478 


0.0054 






2$ lfi3x 


0.2222 


0.0610 


0.0207 








2$ ^ 4x 


0.1111 


0.0198 










Sum 


1.0 


0.64 + 0.05 



Table 2: Probabilities of all possible decay chains for $ pairs in our scenario, "x" in the table 
marks r or, after r decay, r-jet. The reconstructed signatures are marked in bold, the remaining 
signatures were not reconstructed. After omitting the channels that include t —>■ e decay, 69% 
of the total signal is left. In total 64% of the signal has been reconstructed. 

A clear signal extraction from the Standard Model background can be achieved using a set of 
selection rules imposed on a reconstructed event in the following order. 

• SI : events with at least 2 positive and 2 negative muons or jets which have \r]\ < 2.4 and 
Pt > 5 GeV are selected. 

• S2: J2pt (scalar) sum of 2 most energetic positive and negative /i-s or r-jets has to be 
bigger than a certain value (depending on Higgs mass). 

• S3: Z-tagging - if invariant mass of the pair of opposite charged /i-s or r-jets is nearly 
equal to Z mass (85-95 GeV), then the particles are eliminated from the analysis. 

• S"-^ : as $-s are produced in pairs, the reconstructed invariant masses (in one event) have 
to be equal. We have used the condition 

0.8 < mj+/mj- < 1.2. (8) 

If the invariant masses satisfy the condition then we include them to the histogram, 
otherwise we suspect that some /i-s may originate from r decay, and make an attempt to 
find corrections to their momenta according to the method described above. 
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The rule SI is an elementary detector trigger. S2, performing scalar sum of pt, is an untradi- 
tional cut. The advantage compared to the widely used pt cut for a single particle is clearly 
visible from Figure [H The left panel shows that the maximum of Higgs line reaches clearly out 
of the background while on the right panel the maximum is deeply inside the background. 




Figure 1: The left panel shows the distribution of events according to scalar sum of 2 most ener- 
getic (highest Pt) positively and 2 most energetic negatively charged muons or jets (^Pr)- The 
right panel shows the distribution of events considering traditional pt cut for single particles. 
Both figures correspond to luminosity L = 30 fb~^. 



Z-tagging in S3 suppresses pp ZZ and pp — > tiZ background. S4 is based on the equality of 
the invariant masses of like-signed fi-s or r-s. Figure |2] gives a clear picture of the behavior of 
signal and background for the S4 selection rule. Naturally, some freedom is needed due to the 
$ decay width and experimental errors of the detector. We require that the ratio m,f^/m,J~ 
has to be in the region from 0.8 to 1.2. 




Figure 2: Distribution of events according to the ratio of reconstructed invariant masses 
(m$++/m$ — ) (no other cuts are applied). The figure corresponds to luminosity L = 30/6~^. 



While the selection rules SI, S3 and S4 are independent of the Higgs mass, the selection rule S2 
(^Pt cut) has to be optimized for a certain Higgs mass value. The cut may be set to a very 
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high value which ehminates all background events, but inevitable loss in signal may postpone 
the discovery of new physics at LHC Thus it is natural to take the minimal discovery luminosity 
{Lmin) as the optimization criteria. Looking for a cut value that enables to make a discovery 
with the lowest luminosity, we are dealing with small signal and background expectations 
by definition. Simple significance estimators cannot be exploited here. We have used the 
log-likelihood ratio (LLR) statistical method [lU HS] to demand 5a discovery potential to 
be bigger than 95% (1 — CL^+fc > 0.95) as for a discovery criteria. This is a rather strong 
requirement, because it allows to make a discovery (meaning the fluctuation of background 
may mimic the outcome of an experiment with probability less than 2.9 ■ 10"'' (5cr) ) during the 
specified luminosity with a probability of 95% (if s+b hypothesis is correct). The widely used 
convention, that significance should exceed five, gives only 50% discovery potential in Gaussian 
limit and diminishes to very small values when background approaches zero. 

The best value for S2 cut does depend on M$ but is not too sensitive to it. Typically the J2pt 
can be assigned a value with a precision of 100 GeV while affecting the minimum luminosity by 
only a couple of percent. In the Table [3] the approximated middle point of this value is given. 
As the best S2 cut is very strong, it suppresses almost entirely the generated background (being 
combined with the other selection rules). For Higgs masses above 500 GeV the background is 
totally suppressed and the discovery potential criteria meets the requirement for 3 signal events 
(6 invariant masses). Nevertheless we cannot infer that the background is really zero in nature. 
To estimate the statistical error due to final number of generated background events we have 
found 95% upper limit of background according to Poisson statistics (Table HI in brackets). 
Using this limit in LLR analysis we get much higher luminosities for discovery. Even a very 
small background expectancy {b = 0.01) gives some possibility to have one (9.9 ■ 10^^) or two 
(4.9 ■ 10~^) background events in the experiment and these outcomes cannot be interpreted as 
discovery anymore. This phenomenon shifts the minimal required luminocity to much higher 
values denoted as L^^ax in Table El 



Mass of $ (GeV) 


200 


300 


400 


500 


600 


700 


800 


900 


1000 


Optimal for S2 (GeV) 


300 


400 


600 


700 


860 


860 


860 


860 


860 


MC Lrmn (fb-^) 


0.25 


0.93 


2.0 


3.6 


8 


17 


34 


62 


120 


MC Lmax (fb-1) 


0.26 


1.03 


3.1 


7.0 


17 


38 


77 


160 


320 



Table 3: Optimal YIpt cut for different Higgs masses and the corresponding minimal discovery 
luminosities: the lower {Lmin) corresponds to the generated background in our analyses and 
the higher {Lmax) corresponds to 95% upper limit of the background error. 

An example of invariant mass distribution after applying selection rules are shown in Figure [3] 
for M$ = 500 GeV. A tabulated example is given for M$ = 200, 500, 800 GeV in Table S 
corresponding to the luminosity L = 30fb~^. The strength of S2 cut is clearly visible: almost 
no decrease in signal while the number of the background events descends close to final minimum 
value. A peculiar behavior of S4 - reducing the background, while also increasing the signal in 
its peak - is the effect of applying the t fi' correction method described above. 
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I I background S1 
I ] background S2 
background S3 

background S4 
Higgs 500 GeV SI 
Higgs 500 GeV S4 




800 1000 
Inv.mass (GeV) 




800 1000 
Inv.mass (GeV) 



Figure 3: Distribution of invariant masses after applying selection rules (S1-S4) for Higgs 
M$ = 500 GeV and the Standard Model background (L=30 fb"^). The histo gram in the right 
panel is a zoom of the left histogram to illustrate the effects of the selection rules S2-S4. 



Process 


N of invariant masses 




N of $ 


SI 


S2 


S3 


S4 


Energy range 150. ..250 GeV 


M$=200 GeV 


4670 


1534 


1488 


1465 


1539 


tt U 




1222 (168) 


172 (8.5) 


134 (6.9) 


17.6 (3.7) 


tiZ 




21.3 (4.0) 


15.5 (1.0) 


6.3 (1.2) 


2.2 (1.1) 


zz 




95.0 (12.0) 


22.5 (0.7) 


9.8 (0.5) 


1.7 (0.2) 


Energy range 375... 625 GeV 


M$=500 GeV 


119.2 


48.4 


47.5 


46.8 


49.5 


tt U 




178 (28) 


2.1 (0.9) 


1.65 (0.87) 


0.10 (0.35) 


ttz 




6.6 (1.7) 


2.3 (1.0) 


1.0 (1.0) 


0.00 (0.1) 


ZZ 




9.4 (2.9) 


1.4 (0.2) 


0.68 (0.19) 


0.08 (0.09) 


Energy range 600... 1000 GeV 


M$=800 GeV 


11.67 


5.05 


5.00 


4.92 


5.21 


tt U 




77 (12) 


0.00 (0.22) 


0.00 (0.22) 


0.00 (0.07) 


ttz 




2.6 (1.2) 


0.39 (0.4) 


0.39 (0.4) 


0.00 (0.1) 


ZZ 




2.5 (0.8) 


0.34 (0.16) 


0.17 (0.09) 


0.00 (0.02) 



Table 4: Effectiveness of the selection rules for the background and signal. All event numbers in 
the table are normalized for L=30 fb~^. The numbers in brackets mark errors at 95% confidence 
level for Poisson statistics. The signal increases after S4 due to the reconstructed t ^ /i' decays. 



5 Including measurement errors to Monte Carlo 



In this Section we make an attempt to estimate simplified detector effects at the level of Monte 
Carlo analyses. In order to do that we have added overall detection efficiencies for the Monte 
Carlo generated //-s and r-jets - 0.98 and 0.6 respectively. Additionally, we apphed Gaussian 
distortion functions to Monte Carlo produced data for //-s, r-jets and which were used 
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Mass of $ (GeV) 


200 


300 


400 


500 


600 


700 


800 


900 


1000 


Optimal J2pt for S2 (GeV) 


300 


400 


600 


700 


860 


860 


860 


860 


860 


Det Eff Lrmn (fb-^) 


0.526 


1.20 


3.0 


6.6 


15 


30 


60 


111 


200 


Det Eff Lma. (fb-^) 


0.546 


2.19 


6.5 


16.6 


39 


86 


190 


420 


900 



Table 5: Optimal ^ pr cuts and minimal discovery luminosities for different Higgs masses when 
the estimates of detector measurement errors are taken into account. Two boundaries for the 
minimal luminosity are given: the lower {Lmin) corresponds to the generated background and 
the higher {Lmax) corresponds to 95% upper limit of the background error. 

to alter randomly those quantities in the analysis. Although the precision of /i detection is 
sensitive to pt of ^ and \r]\ we use the mean values for a rough estimation. We make the 
following assumptions based on [lH |12[ |13] . 

The direction of muon (r-jet) is altered with the Gaussian distribution: mean /x = 0.0005 
(0.031) and variance cx^ = 0.003 (0.017). The transverse momentum is altered according to 
the PT,rec/PT,MonteCario Gaussiau distribution: mean /i = 1. (0.897) and variance cx^ = 0.03 
(0.089). Both components of missing transverse energy are altered independently according to 
the Gaussian distribution (mean = GeV and variance = 25 GeV) by adding the piece 
to its Monte Carlo value. 

The result of such a distortion is a decrease in both signal and background approximately by 
factor two (Tabled]). As the background and the signal decrease proportionally, the luminosity 
needed for discovery roughly doubles. Remarkably the optimized S2 cut value does not change 
significantly. The proportion of the reconstruction channels in the total analysis has changed 
remarkably as shown in Table El The reason is clearly the small detection efficiency of r-jets. 
The l/u3j channel comprises only 1% of the total signal if the detector effects are considered, 
while in the pure Monte Carlo analysis it forms 9%. The additional possible detector effects 
make the t ^ fi correction even more relevant. 





Percentage 


of the channel after reconstruction 


Decay channel 


4/. 


3/ilj 


2/i2j 


l/i3j 


r — s> /i correction 


Monte Carlo 


21 


28 


26 


9 


16 


MC+ efficiencies 


38 


25 


12 


1 


24 



Table 6: The importance of reconstruction channels at Monte Carlo level and considering 
detector efficiency effects. 



6 Conclusions and outlook 

We have studied possible direct test of little Higgs scenarios which light particle spectrum 
includes a triplet scalar multiplet at LHC experiments. We have investigated the Drell-Yan 
pair production of the doubly charged Higgs boson and its subsequent leptonic decays. In 
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addition to solving the little hierarchy problem, this scenario can also explain the origin of 
non-zero neutrino masses and mixing via the triplet Higgs neutrino mass mechanism. Simple 
connection between the observed neutrino mixing and triplet Yukawa couplings allows us to 
predict the leptonic branching ratios of the triplet. Thus the experimental signatures of the 
model do not depend on the size of the triplet Yukawa couplings allowing direct tests of this 
scenario at LHC 

In our analyses we have considered four /i and/or r final states including up to 3 tau leptons. 
We propose four selection rules to achieve the optimized signal and background ratio. As the 
decays are lepton number violating, we have shown that the background can be practically 
eliminated. In such an unusual situation we have used the LLR statistical method to demand 
5a discovery potential to be bigger than 95% (1 — CL^+b > 0.95) as the discovery criterion. The 
results of optimized cut values are presented in Table [3l Considering the pure Monte Carlo 
study, up to the mass 300 GeV can be discovered in the first year of LHC (L = 1 fb^^) 
and up to the mass 800 GeV can be discovered for the integrated luminosity L = 30 fb~^. 
Including particle reconstruction efficiencies as well as Gaussian distortion functions for the 
particle momentas and missing energy which mimic detector inefficiencies at Monte Carlo level, 
our results show that can be discovered up to the mass 250 GeV in the first year of LHC 
and 700 GeV mass is reachable for the integrated luminosity L = 30 fb~^. 

For further studies of this scenario at LHC progress can be made both physics-wise as well as 
technically. Full simulations of the detector effects are needed which also include the electron, 
muon and tau final states. For better determination of statistical errors coming from the 
background studies bigger SM background datasets must be produced. This requires huge 
computing resources. If these goals can be achieved, the proposed phenomenology opens a 
new window to study the neutrino properties at colliders. In addition to the considerations in 
this paper, one can determine at LHC experiments the hierarchy (normal or inverse) of light 
neutrino mass spectrum, and to estimate the two Majorana phases which are not measurable 
in neutrino oscillation experiments [16] . 
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